Manipulating Topological Edge Spins in One-Dimensional Optical Lattice 
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We propose to observe and manipulate topological edge spins in ID optical lattice based on cur- 
rently available experimental platforms. Coupling the atomic spin states to a laser-induced periodic 
Zeeman field, the lattice system can be driven into a symmetry protected topological (SPT) phase, 
which belongs to the chiral unitary (AIII) class protected by U(l) x Z^ symmetry. In free-fermion 
case the SPT phase is classified by a Z invariant which reduces to Z^ with interactions. The edge 
modes of the SPT phase are spin-polarized, with left and right edge spins polarized to opposite di- 
rections and forming a topological spin-qubit (TSQ). We demonstrate a novel scheme to manipulate 
the spin-polarized zero modes and realize quantum-tunneling driven spin Rabi-oscillations in a con- 
trollable manner. The manipulation of TSQs has potential applications to quantum computation. 

PACS numbers: 71.10Pm, 37.10 Jk, 42. 50. Ex, 71.70.Ej 



Introduction. Since the discovery of the quantum Hall 
effect in two-dimensional electron gas [T], the search for 
nontrivial topological states has become an exciting pur- 
suit in condensed matter physics [2J- The recently ob- 
served time-reversal (TR) invariant topological insula- 
tors (TIs) have opened a new chapter in the study of 
topological phases (TPs), attracting great efforts in both 
theory and experiments [3J 0]. Depending on whether 
the ground states have long-range or short-range en- 
tanglement, the TPs can be classified into intrinsic or 
symmetry-protected topological (SPT) orders [SHZ]- Un- 
like conventional quantum phases characterized by bro- 
ken symmetries, TPs of quantum matter are classified 
by topological orders whose emergence does not break 
symmetry, and states of different topological orders can- 
not change into each other without closing the bulk gap. 
Being protected by the bulk gap, the intrinsic TPs are 
robust against any local perturbations, while the SPT 
phases are robust against those respecting given symme- 
tries [5HH] ■ This property may be applied to the fault- 
tolerant quantum computation [9]. 

While the search of TR invariant TIs has been success- 
ful in solid-state systems, the recent great advancement 
in realizing effective spin-orbit (SO) interaction in cold 
atoms [10H15) opens intriguing new possibilities to probe 
SO effects [16] and TPs in a fully controllable fashion. 
Many theoretical proposals have been introduced in cold 
atoms for the study of TIs [T?rl2"l"] and topological su- 
pcrfluids [22- 25 . Experimental observations of these ex- 
otic phases are, however, delicate due to the demanding 
conditions such as complicated lattice configurations or 
SO interactions. By far the only experimentally realized 
SO interaction [TlTfl"5] is the equal Rashba-Dresselhaus- 
type SO term as theoretically proposed by Liu etal [TU] , 



Therefore, to observe nontrivial topological states with 
currently available experimental platforms is particularly 
desired in the field of cold atoms [2"5] . 

In this letter, we propose to observe and manipulate 
topological edge spins in one-dimensional (ID) optical 
lattice with the SO interaction realizable in recent ex- 
periments [ITULi] . The predicted SPT phase belongs to 
AIII class and is protected by U(l) and chiral symme- 
tries, with spin-polarized zero modes localized on bound- 
aries. Our results not only provide a realistic platform 
to study SPT states in cold atoms, but also have broad 
range of applications including realizing single spin con- 
trol in optical lattice, which is essential for spin-based 
quantum computation. 

Model.— Our model is based on quasi-lD cold fermions 
trapped in an optical lattice, with the internal three- 
level A-type configuration coupled to radiation, as shown 
in Fig. [I The transition \g^) — > |e) is induced by the 
laser field with Rabi-frequency Qi(x) = Qq sm(k x), an d 
the one \g±) — > |e) is driven by another laser f2 2 (a0 = 
f2o cos(kox). I n the presence of a large one-photon detun- 
ing |A| ^> flo and a small two-photon detuning \S\ <C 
for the transitions (Fig. [ija)), the Hamiltonian of the 
light-atom coupling system reads H = H + Hi , with 

ff o = E [^+V a {x)}\g a )(g a \+h8\g i )(g i \ 1 
<t=t,4- 

H x = h^\e){e\-h(n l \e){g^\+n 2 \e){g i \+K.c). 

Here the diagonal potentials Vf^(x) are used to construct 
the ID optical lattice, and o~ yrZ are the Pauli matrices 
in spin space. For |A| » £l , diagonalizing Hi yields 
two nearly degenerate dressed states known as the dark 
and bright states given by [H> \xd) — cos{kox)\g^) — 
sin(A;ox)|5^) and \xb) = sin(fc x)|g t ) +cos(k x)\g l }. Pro- 
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FIG. 1: (Color online) (a) Cold fermions trapped in ID optical 
lattice with internal three-level A-type configuration coupled 
to radiation, (b) Energy spectra with open boundary condi- 
tion in the topological (diamond, F z = 0) and trivial (circle, 
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jecting H onto the subspace of \xb,d) gives rise to the 
equal Rashba-Dresselhaus-type SO interaction which has 
been predicted in theory |10) , and remarkably realized in 
the recent experiments [PD - H5] . In this study we shall 
work in the basis of hyperfine states |<?t,^)- The effect of 
the small two-photon detuning is equivalent to a Zeeman 
field along z axis T z — HS/2, which in experiment can be 
precisely controlled with acoustic-optic modulator. Elim- 
inating the excited state by |e) « ^(f^Lg-f) + ^Ifi^)) 
yields the effective Hamiltonian 



2 

2m 



J2 [V^\x)+T z a z ]\g a )(g a \ 
-=1,2 

[M(x)\g t )( gi \+R.c.], 



(1) 



hVJ. 



where M(x) = Mosin(2fcox) with Mq — ,J ^- represents 
an additional laser-induced Zeeman field. To form a ID 
lattice, we consider the optical dipole trapping potentials 
V^ att {x) = V^ att (x) = -V cos 2 (2k x) with the lattice 
trapping frequency u = (SVok^/m) 1 ^ 2 [55]. This effec- 
tive model describes a system of spin- 1/2 fermions expe- 
riencing a periodic transverse Zeeman field M(x) which 
is commensurate with the ID optical lattice [Fig. [fja)]. 

We next derive the tight-binding model by taking that 
the fermions occupy the lowest s-orbitals <f> sa (a =^,4.). 
From the even parity of <j> S(T , one can verify that the peri- 
odic term M(x) does not couple the intrasite orbitals 4>st 
and <fisi, but leads to a spin- flip hopping given by t l J Q = 

J dx(j)^(x)M(x)(j>[^ (x), as a consequence of the induced 
SO interaction. The spin-conserved hopping is given by 

t, = J dx<f> ( JJ + V]^ a +1) (x). Bearing these re- 

sults in mind we write down the effective Hamiltonian 
in the tight-binding form: H = -t s Y^<i,j>,a Aa^i" + 

E 4 r ^r - %)+[E<i,j> i si4t e J4 + H - c -]' with = 

It can be verified that i«' ±1 = ±(-l)H', 



where tio* = ^? J dxcj) s (x) sm(2kox)(fi s (x — a) with a 
the lattice constant. Redefining the spin-down operator 
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SO j 



C 3i 

H = -U 



e nrxj/a~ we recas ^ the Hamiltonian into 



£ 

<i,3> 



(4 



]T z (h^ - hii) + 



+E 



(0) 



c jt c J-n) 



H.c.]. 



(2) 



We analyze the symmetry of the above Hamiltonian. The 
TR and charge conjugation operators are respectively de- 
fined by T — iKcjy with K the complex conjugation, 

and C : (c ff ,cj.) 1 — > (<7z)<7-<7'(cv> <v). One can check 
that while both T and C are broken in H, the chiral 
symmetry, defined as their product, is respected and 
(CT)H(CT)~ l = H, with (CT) 2 - f. Note the chiral 
symmetry is still preserved if a Zeeman term T y o~ y along 
y axis is included in H. The complete symmetry group 
then reads U(l) x Zj, where U(l) gives particle- number 
conservation and the anti-unitary group Zj is formed 
by {I,CT}- The SPT phase of our free-fermion system 
belongs to the chiral unitary (AIII) class and is character- 
ized by a Z invariant [5HZ] . The H can be rewritten in the 



\d z (k)a z + dy(k)a y ] a)( r>Ck tl 



fc-space H ■■ 

with d y = 2tio sin(fca) and d z = —T z + 2t s cos(fca). This 
Hamiltonian describes a nontrivial topological insulator 
for |r z | < 2t s and otherwise a trivial insulator, with the 
bulk gap E g = min{|2t s - |r,||, 2|t£?|} (Fig. [l] (b)). In 
particular, when T y z — and t s = |iso |j our model gives 
rise to a flat band with nontrivial topology. 

Edge states.— The nontrivial topology can support de- 
generate boundary modes. Considering hard wall bound- 
aries located at x = 0, L, respectively [27] and diag- 
onalizing H in position space H = T~i(xi) with 

H(xi) = -(t s o- z +iUo >u y )c\.c Xi+a +T z a z cl..c Xi +h.c, we 
obtain the edge state localized on left boundary x = as 

M*) = -7Tfi( x +y x ' /a (A-r /a ]ix+>, (3) 

and accordingly the one on a; = L by ipn(xi) — 
^[(A + ) (i - Xl)/a -(A-) (L ~ :E,)/a ]IX-)- Here j\f is the nor- 
malization factor, \x±) — [liil] T are s pi n eigensates of 

a x , and A± = (r„ ± fa - 4t* + 4\t$\ 2 )/(2t s + 2\t£>\). 
Therefore the two edge states are polarized to the oppo- 
site ±1 directions. Note the edge modes ipL,R span the 
complete Hilbert space of one single f/2-spin or spin- 
qubit. Each edge state equals one-half of a single spin, 
similar to the relation between a Majorana fermion and 
a complex fermion in topological superconductors. As a 
result, we expect the robustness of the zero modes to any 
local operations without breaking the U(l) and CT sym- 
metries. These properties of the topological spin qubit 
(TSQ) may be applicable to the fault-tolerant quantum 
computation [9j. Moreover, the Z classification implies 
that single-particle couplings respecting U(l) and CT 
cannot gap out the edge modes in arbitrary ./V-chain sys- 
tem of ID lattices, while, interestingly, we have confirmed 
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that weak interactions can break Z down to Z4 classifi- 
cation 28 . This result suggests an interesting platform 
to study classification of SPT phases with interactions. 
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FIG. 2: (Color online) (a) Wave functions for zero modes 
\iPl,r)', (b) 1/2-particle fractionalization (seen by AJVi^) for 
zero modes. The parameters tao = 0.4t s and F z = 0.3t a , with 
which the localization length of bound modes £o=2.36a. 

Existence of zero modes leads to particle fractionaliza- 
tion, which is another direct observable in experiment. 
A zero mode is contributed half from the valence band 
and half from the conduction band. Therefore, an edge 
state carries +1/2 (—1/2) particle if it is occupied (unoc- 
cupied) |28j . This result is confirmed by numerical sim- 
ulation shown in Fig. [2j where we calculate ANi^{xi) = 
Jq % dx' Ni >2 (x') — Xi/a at half-filling, with Ni( 2 )(x) the 
density of fermions when the left (right) edge mode is 
filled. The fermion number carried by an occupied (un- 
occupied) edge mode is then given by ni( 2 ) = AiVi(2)(£) 
with £ £o- Here £0 = — a/ln|A+| is the localization 
length of ipz R- The 1/2-fractionalization is clearly seen 
when £ is several times greater than £0 (Fig.[2](b)). Being 
a topological invariant, the 1/2-fractionalization can be 
confirmed to be robust against weak disorder scatterings 
without breaking the given symmetries. 

Correlation effects.— A particular advantage in cold 
atoms is that one can investigate correlation effects on the 
predicted SPT phase by precisely controlling the interac- 
tion. For spin- 1/2 cold fermions the onsite Hubbard in- 
teraction UJ2i n it n il can be well controlled by Feshbach 
resonance [26]. We adopt here the Abelian bosonization 
approach combined with renormalization group (RG) 
analysis to probe interaction effects in single-chain sys- 
tem at critical point. With low energy approximation 
we linearize the fermion dispersion near Fermi points 
k F = Ifc by c s (x) « X)r,a=±i ip rs e lrkFX , where r/s is the 
chiral/spin index. Using the standard bosonization for- 
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e V2 



[r</> p -0 p +s(r<l>„-6„)] 



we reach 



the bosonized Hamiltonian density Ti 
withH, = ^[K^d^f+K-^d^ 2 



) 2 H 



Yl lv =p, a ^-\-'rL\ 
^cos2V2>„ 



sin y/2(f) p cos \JlQ a 



(27T0) 

cos y/2(j)p sin y/28, 



and 7ii = 

Here g p><7 arise from umklapp scattering and spin 
backscattering, respectively. The bare values are v v — 



2t 



(0) 



w = Ty in units of t s 



At 



v F K Vl g p = -g a = U,u 

half-filling Ti respects 0(4) S SU(2) x SU(2) symmetry if 
u = w = [3U]. For U > 0, the umklapp scattering term 
cos 2\^2(f>p drives the system into a Mott insulating phase 



with a charge gap. The spin sector remains at the fixed 
point with K* = 1, described by SU(2)i Wess-Zumino- 
Witten model 3 lj . When u,w 7^ 0, such symmetry is 
reduced to the U(l) subgroup. 

For the large-?/ limit, the charge boson <f) p is pinned 
to its classical minima cos2v / 2</> P = 1, and the system 
is a Mott insulator, regardless of the mass terms. In the 
weak-coupling regime, the interaction effects on the phase 
boundary can be studied with RG method up to one- 
loop order [2J. The RG flow equations for the coupling 
constants u,v,g p ,g& are given by |28) : 



dit 
~dl 
dw 
~dl 

dl 



K n 



3-K n 



9 P u 
9 P w 



9aU 
grrW 



9l 



dff CT 
dl 



Airvp Attvf 
9l 



(4) 



Here I is the logarithm of the length scale. The renormal- 
ization of Luttinger parameters K v has been neglected as 
it is a higher order correction. For U > 0, g a marginally 
flows to zero so we drop it off below. g p is marginally 
relevant and can be solved by g p (V) — ^^^w . Sub- 
stituting this result into RG equations of u and w yields 
after integration u(l) = u(0)[l-2^]*e( 3 - K >-W 2 ,w(l) = 

w(0)[l + ^ 1 }^e {3 - K " )l/2 . The physics is clear: the re- 
pulsive interaction (g p > 0) suppresses SO induced mass 
term u while enhances the trivial mass term w. The fate 
of the system depends on which of u and w reaches the 
strong-coupling regime first. Assuming |g p (0)^| <C vp, we 
find the TP transition occurs at 



u(0)=[w(Q)]\ 7 «1 



9p(0) 



4tt?; f (3- Kp) 



(5) 



This gives the scaling law at the phase boundary with 
interaction. Note 7 < 1 for U > 0. The above scaling re- 
lation implies that a repulsive interaction suppresses the 
SPT phase. Accordingly, if initially the noninteracting 
system is topologically nontrivial with m(0) > w(0) > 0, 
increasing U to the regime w(0) < [u>(0)] 7 drives the sys- 
tem into a trivial phase. 

Single spin control. —Now we study an interesting ap- 
plication of the present results to realizing single spin 
control. Besides the edge modes localized on the ends, 
TSQs can also be obtained in the middle areas by cre- 
ating mass domains in the lattice. This can be achieved 
by applying a local Zeeman term T y or T z . For exam- 
ple, we consider T z — everywhere, but r„ = r for 
Xi < x < X2 and Y y = otherwise. The local T y can 
be generated by applying another two lasers which cross 
with the ID lattice and couple the atoms in the area 
xi < x < X2 to induce a local resonant Raman coupling 
between \g^) and \g±) (Fig.[3ja)). Employing a 7r/2-phase 
offset in the Rabi-frequencies of the two lasers, the Ra- 
man coupling takes the form Toa y , with Tq controlled 
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by the laser strength. When \T \ > 2\tso \ a mass do- 
main is created, associated with two midgap spin states 
\tjj±) respectively localized around x — X\p (Fig. |3ja)). 
The width Ax = x 2 — x\ and height of the domain are 
respectively adjusted by the waist size and strength of 
the two laser beams. Due to the nonlocality of the TSQ, 
creation of a single qubit here is not restricted by the 
size of the laser beams. This is a fundamental differ- 
ence from creating conventional single qubit by optical 
dipole trapping which requires tiny-sized laser beams to 
reach a very small trapping volume |33) . Note in realistic 
case the laser induced T y may vary fast but not in the 
form of step functions around x = Xx,%, which, however, 
does not affect the main results presented here. Cou- 
pling between and \ip-} results in an energy split- 
ting 28 oc e -(l r °l- 2 l*-l)Az/(2at s ) ; w hich is controlled by 
F and Ax (Fig. |3ja), lower panel). In the limit (|r | — 
2|£so \)Ax/(2at s ) » 1, such coupling is negligible, and 
the two zero modes consist of a single spin qubit which is 
topologically stable. Let be initially occupied while 
\ip~) be left vacancy. Reducing | To | smoothly can open 
the coupling in \ip±) and lead to spin state evolving as 

m MO) = <*(t)<P-(x - *i)IX-) + 0(t)<P+(x ~ x 2 )\x+), 
with a(0) = 0, (5(0) = 1, and ip± the spatial parts of the 

bound state wave-functions. The spin-polarization den- 
sities are given by s x , y , z (x,t) = (ip(x)\a Xt y tX \ip(x)) , and 
the spin expectation values S x>VtZ (t) — J dxs Xi y tZ (x,t). 
It can be verified that S y (t) = 0, and 

S x (t) = \a(t)\ 2 -\m\ 2 , 

r i (6) 

S z (t) = 2Re[a{t)/3*{t) / dxip+(x)tp_(x)]. 

This phenomenon is analogous to spin precession with the 
rotating angle yielding j(t) — 2 f Q dt'8(t'). We have then 
a = cos -f(t) and j3 = sin 7(f)- The amplitude of S z (t) 
is given by S™ a * = \ J dxtp*,(x)ip-(x)\, which measures 
the overlapping integral of tp±. Accordingly, if we ap- 
ply the local Zeeman field along z rather than y axis, we 
shall obtain that the spin evolves in the x-y plane. Note 
the spin Rabi-oscillation is induced by quantum tunnel- 
ing. Therefore it is associated with a tunneling current 
given by J m (t) — ~^ff d t \a(t)\ 2 between x\ and x 2 - In 
experiment the internal states of a single atom can be de- 
tected without energy transfer [35], which is applicable 
to observe the spin Rabi-oscillations, while the oscillation 
of S x (t) can be more conveniently observed by measur- 
ing the number of fermions (n±(t)) localized on x± t 2, and 
Jm(t) can be detected by measuring the change rate with 
time of such fermion numbers. 

We show in Fig. [3] (b-d) the numerical simulation 
for single spin control with the parameter regime that 
t s = 3.15kHz, t$ = 0.4t s , and Ax = 10a. For t < 0, 
r = 8is? and the coupling in \ip±) is negligible. Re- 
ducing r at t > leads to spin evolution and by fixing 
Tq — 2.5iso for t > t\ the spin oscillates with a period of 




t (ms) * (ms) 

FIG. 3: (Color online) Spin Rabi-oscillations with the param- 
eters t s — 3.15kHz, tlo = 0-4:t s , and Aa; = 10a ~ 4fim. (a) 
Mass domain created by setting |Fo| > 2ti'o l for xi < x < X2 
which localizes a spin-qubit composed of two bound modes 
V-'i) on x = 21,22, respectively; (b) Spin Rabi-oscillation by 
smoothly reducing \Tq\ from 8iao to 2.5ti® ; (c) The mass cur- 
rent J m (t) and expectation values of particle numbers (n±(t)) 
in states \tp±)\ (d) Spin-flip operation by controlling that 
t(*s) ~ (2 m + l) 71 " with m — 1. The initial spin state |x+) 
(Points A) flips to be |x-> (Points B). 



5.984ms (b-c). Note the quantum state of the spin can be 
precisely controlled by properly manipulating j(t). For 
example, in Fig. [3] (d) we demonstrate the spin-flip op- 
eration \x+) ~^ \X—) by requiring 7^3) = (2m + 1)tt. 
Here m £ Z and in (d) we take m = 1. Note one may 
integrate multiple TSQs with e.g. atom-chip technology 
and individually control them by creating multiple mass 
domains in the ID lattice. The precise manipulation of 
such integrated TSQs may have interesting applications 
in developing scalable spin-based quantum computers. 

Before conclusion we estimate the parameter values 
for realistic experimental observations. For example, in 
40 K atoms we have the recoil energy En/h = hk 2 /2m = 
48kHz using red-detuned lasers of wavelength 780nm to 
form the optical lattice. Taking that Vq = 5Er and 
M = 2Eji, we have that the lattice trapping frequency 
lu = 214kHz, and hopping coefficients t s /H ~ 3.15kHz 
and /h ~ 1.3kHz. Then the bulk gap equals E g /h — 
2.6kHz for T z = 0, indicating a temperature T = 19nK 
for the experimental observation. 

Conclusions.— hi conclusion, we have proposed to ob- 
serve and manipulate the topological edge spins in ID 
optical lattice with currently available experimental plat- 
forms. In the free-fermion case the SPT phase belongs 
to the AIII class characterized by a Z invariant which 
breaks down to Z^ with interactions. The spin-polarized 
zero edge modes in the nontrivial SPT phase form TSQs. 
Our results are of broad range of interests and applica- 
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tions including the manipulation of TSQs which may be 
applicable in developing spin-based quantum computers. 
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Manipulating Topological Edge Spins in One-Dimensional Optical Lattice — Supplementary Material 

In this supplementary material we provide the details of some results in the main text. 



Appendix A: Topological Classification 
1. Noninteracting regime 

With the inclusion of both T y and T z , the generic Hamiltonian obtained in the main text is given by 

H = -£[*«(c5t^'+ 1 t ~ S'A+u) - 4o(CjfC j+n - ctfCj-u) + H.c] + ^[r z (n jt - - T y (ic] tCjl + H.c.)] 

3 3 

= EC^^V' ( A1 ) 

fc,(T 

where "H(fc) = d y (k)a y + d z (k)a z with d y = T y + 2t$ sin(fca) and d z = T z — 2t s cos(fca). The full symmetry of 
the system is U(l) x (CT), where the phase transformation operator U(l), charge conjugation operator C and time 
reversal operator T are defined as: 

U(6)cU(e)- 1 = e i6 c, CcC' 1 = aj; 
TcT~ x = -ia y Kc, (CT)c(CT)- 1 = o x &. 

The following commutation relations can be checked: C 2 = -T 2 = (CT) 2 = 1, {C,7} = 0, and [CT, C/(l)] = 0. The 
subgroup {I,CT} is anti- unitary and can be denoted as Z% , and the complete symmetry group can also be written 
as U(\) x Z% ■ Owning to this symmetry group, the free fermion system belongs to chiral unitary (AIII) class and is 
characterized by a Z invariant in the noninteracting case. 

If the Zeeman field along y axis vanishes, i.e. T y = 0, alternatively the symmetry group can be chosen as U(l) x 
x(C x T) where C/(6')cC/(0)~ 1 = e t9 c, CcQ- 1 = a x c\ Tc^ 1 = Kc, with C 2 = V 2 = (CT) 2 = 1, and [C,T] = 
1,777(0) = U(—8)T,CU(9) = U(—6)C. In this case both T and C are symmetries of the Hamiltonian and the system 
then belongs to the BDI class which is also classified by a Z invariant in the non-interacting case. 



2. Interacting regime 

For a N-chain system under the half-filling condition, the total ground-state degeneracy is C^ N without interactions. 
Let c m . L / R and $ mL i R be the annihilation and creation operators of the left/right edge mode for the ra-th chain, 
respectively, where cl = (cf + cj.)/\/2 and cr = (cf — c±)/y/2. The edge states of the m-th chain can be written as 
c' m L \0) m and R \0) m , where |0) m is the ground state for the bulk. 

An interesting question is what happens if we turn on interactions. It turns out that in the presence of interactions, 
the Z classification breaks down to Z 4 . To confirm this result, we will study the ground states of a N-chain system 
step by step. 

First, for N — 2, a generic U(l) x (CT)-symmetric interaction between the edge zero modes reads 

H12 = [ F 12 4) Cl, s Cl, S 4, s C2, s + (Vi2 4) )*Cl, S ct )S C2, a 4 ia 

s=LM 

= ^Wi.aAAs - - 4, a. + !)• ( A2 ) 

s 

With above interaction the minimum degeneracy of the 2-chain system is two-fold, which is obtained when V$ = 
(y$)* = -V < 0. It is straightforward to check that the states |01)i|01) 2 and |10)i|10) 2 have energy -2V while 
the states |10)i|01)2, |01)i|10) 2 ,|00)i|ll) 2 , and |ll)i|00) 2 have a higher energy 0. Then the ground state is two-fold 
degenerate, as given by |01)i|01) 2 and 1 10) 1 1 10) 2- 
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Secon d, for N = 3, the ground state degeneracy Cf = 20 can be reduced to 2 by two-chain interactions according 



to (|A2j): H ult = H\2 + i?23- For instance, if V$ = Vjs = —Vq < 0, the two-fold degenerate ground states are 



and if Vu 



|01)i|01) 2 |01) 3 , |10)i|10) 2 |10) 3 ; 
-V23 = — Vq < 0, the ground states read 

|01}i|01}a|10> 3 , |10)i|10) 3 |01) 3 . 



(A3) 



(A4) 



However, for the 3-chain system, one should also consider three-body interactions. One of the possible 3-chain 
interactions reads 



#123 = X!( y i23 a lAs4,As4,As + ^123Cl,*4,A«4,A«4,J« 



(A5) 



where V^l = (^123)* • ^ is straightforward to see that above interaction is identical to a summation of two-body 
interactions and can not split the degeneracy of the ground states. Other possible 3-chain interactions include 



#123 = 



4 s C 2 s C 2,sC3,s + Cl, s C2, s C. 



s 

^(2cU-l)(4/3, 



2,s c 3,s 



h.c. 



u 3,s 



c l,s), 



(A6) 
(A7) 



It can be easily checked that in the case ( A3 ) the perturbation H' 12 3 have zero matrix elements in the ground-state 



subspace and hence can not spilt the degeneracy. The condition is similar in another case (A4). The difference is that, 



in case (A4|, H' 1Ti mixes the ground states with higher energy states. For instance, the state |01)i|01)2|10)3 is mixed 
with |ll)i|01) 2 |00) 3 to lower its energy. At the same time the state |10)i|10) 2 |01) 3 is mixed with |00)i|10) 2 |ll) 3 to 
lower the energy with the same amount. As a result, the two new states, as the new ground states, are still degenerate. 
Therefore, the topological properties of 3 chains are stable against interactions respecting the symmetry, and we need 
to investigate the 4-chain system. 

Finally, for N = 4, it turns out that we can find a path to smoothly reduce the ground state degeneracy to 1. It is easy 
to verify that degeneracy of the ground states can be reduced to 4 under two-body interaction _ffi nt = H12 + -ff 23 + H34 
with V12 = —V23 = V 34 = — Vq < 0. These four ground states are given by 



|01)i|01) 2 |10) 3 |10) 4 , |10) 1 |10) 2 |01) 3 |01) 4 , 
|ll}i|ll}2|00> 3 |00>4, |00) 1 |00) 2 |11) 3 |11) 4 , 



(A8) 



with the energy equal to — 4Vo- However, the above ground states can be further gapped out by taking into account 
the following interactions: 



#1234 = 53^1234 C L4, S C 3, S C 4 , s + V^CM C 2 , s c\ s c\ s ). 



(A9) 



In the 4-dimensional Hilbert space spanned by the ground states in Eq. ( A8 ) , the above interaction can be written in 
the matrix form 



#1 



234 



(4) 



1234 



fo 1 l\ 
11 
110 

\1 1 0/ 



(A10) 



whose eigenvalues are —2^234,0, 0, and 2^234- We therefore obtain the single non-degenerate ground state with the 
energy —2| 1^234 1 —AVq. This implies that under interaction the 4-chain system can be smoothly connected to a trivial 
phase without closing the bulk gap, and we therefore complete the proof that the Z classification can be broken down 
to Z 4 under interactions. 



Appendix B: Particle fractionalization 

We prove in this section that each edge state leads to 1/2-fractionalization. A convenient way is to consider the 
semi- infinite geometry which has the open boundary at x — 0. We then calculate the particle number of the zero 
mode localized on this boundary. Note the total number of quantum states in the system is given by 

N= ^2{tp E \h E \4>E) + ^2{ip E \h E \^E) + {MfiEHa), (Bl) 

E <0 E >0 

where wc denote by hp = I the state number operator and \ip E ) is the eigenstate with energy E. Since the Hamiltonian 
satisfies {H,a x } — 0, the energy spectrum is symmetric. We have then ^2E<o^E\n E \ip E ) — X)_b>o(' ! /'£;|^-b|' ! /'b)- It 
follows that 

Y,^E\n E \^E) = \[N- {Mn E \^o)]. (B2) 
E <o 

The particle number of the zero mode depends on its occupation. If the zero mode is unoccupied, the particle number 
of it is given by 

n = ^ [(V'-e|"-b|V'-b)i - (ip E \n E \ip E )o] . (B3) 
E <o 

Here ()i and ()o represents the cases with one (topological phase) and zero (trivial phase) bound modes, respectively. 
Using the Eq. ( B2 ) one finds directly 

no = -^(ipo\n E \^o) = ~\- (B4) 

Similarly, if the zero mode is occupied, the particle number is no ~ 1/2. It is trivial to know that this result can 
be applied to the case with two boundaries located far away from each other, say respectively at x — and x = L. 
Since the two zero modes are obtained independently, each of them carries +1/2 (—1/2) particle if it is occupied 
(unoccupied) . 



Appendix C: Derivation of the RG equations 



In this section we provide the derivation of the one-loop RG equations. We consider that in Hamiltonian ( Al ) only 
one Zeeman term, e.g. T y is nonzero (The case with T z ^ can be studied in the similar way). We find it convenient 
to redefine that Cjj, — > e mx il a Cj\ r and rewrite the Hamiltonian ( Al ) in the form: H = —t s Yl<i j> a ^ta^jc + -^SO +-Hzj 

where H so = ti°' £V ;{-l) 3 (c^c J+ i^ - c^Cj-i^ + h.c.) and H z = T y ( ic ]i c 3t ~ n - c -)- The low-energy physics 

can be well captured by the continuum approximation (x — ja): 

c ja « MipRAx)e lkFX + ipL*(x)e~ lk n, (CI) 

with kp = 7r/2. The continuum representation of the two mass terms are given by: 

H so wiu fdx {^ L a x i} R - ^ R o x i> L ), H z =w [dx {^ R a y ^ L + ^ L a^ R ). (C2) 



Here u — 2t SOl w = T y . Using the standard bosonization formula ip rs = e ^ r ^ p e p+ s ("^ w jt n r — R,L 
and s =f, i, we reach that 



"Hso = — sin V2(/) p cos \f26 a , H z = — cos V2</>„ sin \p28 a . (C3) 
Tra 7ra 



The allowed four-fermion interactions in the (unperturbed) Hubbard model at half-filling is highly constrained by 
the §U(2) x SU(2) symmetry. It is convinient to define current operators: 

Jr = V^aVVa, Jr = ~Z Ca^r^ h = VVa Cot/3 (C4) 
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Here e is the fully anti-symmetric tensor. The general form of four-fermion interaction is given by 

H int = 2gJ R ■ J L + ^(I^Il + h.c). (C5) 

We then derive the RG flow equations to understand the fate of the topological phase transition driven by the 
competition of the staggered Zeeman term and the spin-orbit coupling term. The tree level term can be easily read 
off from the scaling dimensions of Hso and Hz in their bosonization form ( |C3| , both of which are (1 + K p )/2. Thus 
it is necessary to go to the next order in the perturbative expansion. The one-loop order RG flow can be most easily 
derived by calculating the operator algebra of the various operators [U |2] . 

To derive the RG equation, we consider the partition function in the Euclidean functional integral representation 



Z = J VipVipe-*. (C6) 

The Euclidean action S = J dxdr (Hq+Hi) where Ho is the unperturbed Gaussian part and the Hi — Hso+Hz is the 
perturbation. To perform the RG, we expand the exponential to the second order in H\. Let us write H\ = J2giOi. 
Then the second order term is given by 



(C7) 



We have introduced complex coordinates z, w where z = vt — ix. Here we assume the following operator product 
expansion (OPE) 

0^)0,(0) ~ -r^^Ok + regular terms, (C8) 

which is sufficient for our purpose. The OPEs are valid when two points z and are brought close together, as 
replacement within correlation functions. 

At this point the cutoff prescription needs to be carefully specified. We will choose a short-distance cutoff a in 
space, but none in imaginary time. For a rescaling factor &, we must then perform the integral 

/OO y 27T 
dT^— 2 ^ = —lnb. (C9) 



Eq. (C7) then becomes 



'-g igj lnb f O k . (CIO) 
Upon re- exponentiating we obtain the one-loop RG equation 



47TV* 



^ = -r^Ew^- ( cn ) 

Here I = In b. 

Now let us be more specifit. Besides the current operators defined in ([C]), we also need to define 

Mi =^Vl, i = x,y. (C12) 
The fermionic field operators satisfy the following OPEs: 

^(*)^(0) ~ ^, «*)4,(0)~^. (C13) 

The OPEs of the currents and the mass Mj can be calculated from ( C13[ ) using Wick's theorem. Those between the 
currents are standard and can be found in [2] which we do not duplicate. Below are the needed ones: 

Mi(z)J R ■ J L (0) ~ — -^M t (z), Mj(z)J R ■ J L (0) - jAr^M}(z) 
16tt 2 |z| 2 16tt 2 |z| 2 

Ml (z)[I R ll +h.c.}(0) - -^^Mjiz), Mt( Z )[l4+h.c.](0)~~M 1 ( Z ). 
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Applying the above formalism to the model at hand, we find 



du _ 3 - K p _ g p u g a u 

-11 — ~r 



dl 2 Attvf ^7tvf' 

dl 2 4ttv f 4ttv f ' v ' 

dg<j = _J_ 2 % = _J_ 2 

dl ttvf al dl ttvf 9 



If we take the interaction in the lattice model to be an on-site Hubbard form: Hjj = Uj^i n i^ n iii the bare values of 
the coupling coefficients are 

9 P = U,g a = -U. (C16) 

Then for the repulsive interaction U > 0, the coupling parameter g a of the spin sector is marginally flows to zero, 
while g p is marginally relevant. 
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